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SUMMARY 


Important scientific questions concerning pure fluids and fluid 
mixtures near critical points are identified and are related to the 
progress of several disciplines. Consideration is given to questions 
about thermodynamic properties, transport properties, and the complex 
nonlinear phenomena which occur when fluids undergo phase transitions in 
the critical region. We discuss, quantitatively, the limitations to 
answering these questions by experiments in the earth's gravitational 
field. The distinction is made between practical limits which may be 
extended by advances in technology and intrinsic ones which arise from 
the modification of fluid properties by the earth's gravitational field. 
The kinds of experiments near critical points which could best exploit 
the low gravity environment of an orbiting laboratory are identified. 
These include studies of the index of refraction, constant volume 
specific heat, and phase separation. 

Key words: Critical point; gravity effects; phase transitions. 
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1. Phase Transition and Critical Phenomena in Fluids : Scientific Questions 

and Interactions with Other Disciplines . 

1.1 Introduction 

The justification for space experiments in the area of "Phase Transition 
and Critical Phenomena in Fluids" is essentially scientific. This particular 
area of science has strong interactions with solid 9 tate physics, fluid physics 
and other areas of physics. We may expect that scientific advances in these 
areas will bring, indirectly, long terra technological and economic benefits; 
however, ^t cannot be said that the space experiments in theraselves will produce 
such benefits. Therefore the justification for space experiments in 
phase transition and critical phenomena in fluids depends upon first establishing 
that this is an important, active field of scientific research and secondly, 
upon establishing that space experiments are likely to have a major impact on 
its further development. A large part of the Importance of this area of 
research results from its initimate relationship to many areas of solid state 
science. This will be discussed in the next section. In order to develop a 
rationale for the argument that space experiments can have a major impact in 
this area of research, we first discuss the important unanswered questions on 
the subject; we then describe how gravity imposes severe limitations in 
finding experimental answers to these questions; and, finally, we indicate to 
what extent these limitations can be reduced by conducting experiments in 
space. 
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1.2 The Relationship of "Phase Transition and Critical Phenomena In 

Fluids'* to other Disciplines . 

Many of the research subjects described as Included In research on 
"phase transition and critical phenomena In fluids" are Intimately related to 
similar subjects In solid state physics. The experimental phenomena observed 
and the theoretical Ideas used to describe both fluids and solids have a high 
degree of overlap. Indeed, this overlap is recognised by the "Phys.cs and 
Astronomy Classification Scheme. 1975" adopted by the Abstracting Board of the 
I teroatlonal Council of Scientific Unions. This classification scheme, 
which is used for grouping scientific articles for publication and Information 
retrieval, groups "phase equilibria, phase transitions and critical points" 
together for both solid and fluid systems. We expect that important experiments on 
fluid systems, such as those to be defined in this study, will have a substan- 
tial impact on corresponding areas of solid state science. The basis for our 
expectation is the strong interrelationship in recent history between the 
study of phase transition and critical phenomena In fluid and in solid systems. 

The interrelationship arises from enalogies between many phenomena which occur 
near the critical point of the liquid-gas phase transition In pure fluids and 
phenomena which occur at a corresponding "critical point" of phase transitions 
in very different kinds of physical systems. A list of these different kinds 
of physical systems would Include ferromagnets near the Curie point, anti- 
ferromagnets near the Ndel point, certain ferroelectrlcs near the ferro- 
electrlc-paraelectrlc phase transition, alloys exhibiting order-disorder 
transitions (e.g. 6-brass), crystals exhibiting order-disorder structural 
transitions (e.g. ammonium chloride), binary liquid mixtures near consolute 
points, multicomponent fluid mixtures near plait points, and helium near its 
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superfluid transition. (For general references on these subjects see: 

Elcock, 1956; Stanley, 1971; Heller, 1967; Fisher, 1967). 

In each of these physical systems, there is a "quality" (such as the 
difference between the densities of the liquid and vapor phases of. a fluid) 
which distinguishes between two phases which coexist in thermodynamic equilib- 
rium. As the critical point is approached this quality gradually disappears. 
The gradual disappearance of this quality is the most obvious feature which 
distinguishes phase transitions with critical points from other sorts of phase 
transitions such as melting, sublimation, or crystallographic phase changes. 

The diverse physical systems with critical point9 have a number of 
experimental properties in common near their respective critical points. They 
all exhibit a specific heat anomaly, long thermal relaxation times, a marked 
sensitivity to external fields and impurities, and the quality which is 
disappearing shows large long-lived fluctuations. The density fluctuations in 
a pure fluid are responsible for the striking visual effects called ciitical 
opalescence. Analogous fluctuations in solid systems manifest themselves in 
ways such as enhanced neutron scattering, ultrasonic attenuation, and 
electrical resistivity. 

The existence of diverse physical systems displaying analogous critical 
phenomena provides the experimenter with the opportunity to choose both the 
best material and the best technique for studying one or another aspect of 
critical phenomena. For example, the size of fluctuations is rather easily 
studied in magnetic solids with neutron scattering, even very far from the 
critical point. The same technique cannot be used extremely close to the 
critical point because the fluctuations become much larger than the wavelength 
of the neutrons available. On the other hand, fluctuations in fluids may be 


4 






studied easily quite close to the critical point using scattered light. We 
wish to emphasize that extremely close (often near 99.9% of the 
critical temperature) to critical points, all experiments on solids are limited 
by lattice strains which result from impurities, vacancies, etc. Thus 
experiments with fluids (which continuously "anneal" themselves) offer the best 
opportunities for closest approach to the critical point. 

Since critical phenomena occur in a wide range of systems, one might expect 
that their theoretical explanation does not depend upon the detailed nature of 
the Interatomic forces in each system but rather, c \>ld be based upon an 
accurate treatment of interactions which contain only those few general 
properties of the true interatomic forces which are needed to make critical 
points occur. Accordingly, similar theoretical treatments of microscopic 
interactions have been used to interpret data on these diverse systems. The 
van der Waals model of a fluid, molecular field model of a ferromagnet, and 
the Bragg'-Williams model of a binary alloy are closely related in approach and 
yield similar predictions concerning thermodynamic properties near the critical 
point of each of the corresponding systems. There is an exact correspondence 
between the lattice gas model of a fluid and the Islng model of a ferromagnet. 
This correspondence has been exploited as a plausibility argument for taking 
a wide variety of theoretical results based on lattice models of critical 
phenomena and applying them to describe real fluids as well as real solids 
near critical points. Interestingly, several important recent advances in the 
theory of lattice models are outgrowths of the "renormalization group" 
technique, a technique developed for problems in theoretical high energy 
physics (Ma, 1973). 
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Other aspects of phase transition phenomena In liquids have solid state 
counterparts. Away from the critical point Itself, phase transitions in both 
fluids and solids may be initiated by nucleatlon processes or may occur 
spontaneously through a splnodal decomposition mechanism. The theoretical 
description of these processes make use of the concepts of interfacial energy, 
bulk free energy, and diffusion. The concepts are applicable to both fluid 
and solid systems. The time scales for phase changes are very different for 
fluid and solid systems, thus facilitating complementary experimental studies 
which benefit the understanding of both fluids and solids. (Cahn, 1968; 

Langer and Bar-on, 1973; Schwartz et al., 1975). 

In summary, basic scientific studies of phase transitions and critical 
phenomena in fluids are closely related to the study of similar phenomena in 
solid systems. It is also quite likely that future developments in chemical 
engineering will exploit the progress now being made in the understanding of 
thermodynamic and transport properties of fluids near critical points. In the 
study of pure fluids, equations of state for the critical region have been 
developed recently which are much more accurate representations of data than 
are standard engineering equations. (Levelt Sengers et al. 1976). These new 
representations have very few parameters which must be adjusted for each fluid 
thus, they require fewer experimental measurements for reliable predictions. 
The extension of these new equations to special mixtures has been demonstrated 
and work on extending them to mixtures of engineering Interest is in progress. 
(Leung and Griffiths, 1973; D'Arrlgo et al. , 1975). Similar .advances have 
occurred and are occurring in the correlation and prediction of transport 
properties in the critical region. Thus we can expect that advances in the 
scientific understanding of phase transition and critical phenomena in fluids 
will influence chemical engineering practice. 
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1.3 important Questions in Phase Transition and Critical Phenomena in 

Fluids . 

We will consider some of the Important scientific questions in the 
area of phase transition and critical phenomena in fluids which are amenable 
at least to a partial answer by experiment. Naturally this cannot be done 
without some reference to current theoretical ideas. We will first discuss 
pure fluids and then fluid mixtures. Ferhaps the single most important 
question in this area is: exactly how much alike are the superficially 

analogous phenomena occurring in the large variety of systems showing critical 
points? Rigorous renormalization group calculations which apply to large 
classes of model systems indicate that static correlation functions and all 
the thermodynamic properties (each of which can of course be calculated from 
the correlation functions) are "universal” (in the sense that the same 
description applies to each) asymptotically close to the critical point. More 
specifically, they predict that the correlation functions will depend upon the 
dimensionality of the system under consideration (e.g. whether it is two 
dimensional like a membrane or three dimensional like a crystal) and the 
dimensionality of the "order parameter" (or the "quality" which vanishes as 
phases become identical at the critical point). Thus the correlation 
functions for a fluid (where the order parameter is the scalar difference in 
density between coexisting phases) will differ from correlation functions for 
isotropic magnets (where the order parameter is a three-dimensional vector) . 
(Wilson and Kogut, 1974). 

Variables such as the lattice structure (e.g. hexagonal close-packed or 
cubic) and the presence or absence of second nearest neighbor interactions 
etc. are expected to be Irrelevant in determining the functional 
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form of the asymptotic expansions of correlation functions and thermodynamic 
functions; however these other variables will determine the numerical value 
the critical temperature itself and numerical values of various amplitudes in 
the expansions for correlation and thermodynamic functions. At the present 
time, nearly all thermodynamic experiments on pure simple fluids indicate small, 
but experimentally significant, differences between the measured properties and 
those calculated for three dimensional lattice models with a scalar order 
parameter. (Levelt Sengers and Sengers, 1977). It is possible that these 
differences indicate that fluids do not belong to the same "universality" class 
as the lattice models with a scalar order parameter. It is also possible these 
differences indicate that the amplitudes of the corrections to asymptotic behavior 
are quite different for fluids than for the lattice models studied. An experimental 
distinction between these two possibilities would be of great value in understanding 
the range of applicability of this important idea of universality. It is possible 
that this distinction could be made by measurements of thermodynamic properties 
closer to the critical point than is now possible in experiments carried out 
in the earth's gravitational field. The v< y same situation exists to a lesser 
degree when different pure fluids are intercompared. The apparent differences 
between fluids are smaller (hence subject to greater experimental uncertainties) 
than the differences between fluids as a class and lattice models; however, the 
same question exists regarding their origin. Again, the answer could result 
from experiments carried out closer to the critical point than presently 
possible. 

The decay of the range of density fluctuations (or the pair correlation 
length) may be measured by measuring the angular distribution of light. 

X-rays, or neutrons scattered by a fluid. The temperature and density 
dependence of this quantity is predicted to be universal, but a significant 
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test of universality cannot be made until tv»G conditions are met: 1) The 

wavelength of the incident radiation must be short compared with the range of 
correlation. 2) The wavelength of the incident radiation must be long '-red 
with the range of the interatomic forces. Neutron and X-ray scatter s 
experiments easily satisfy the first criterion; however, the second criterion 
is not well satisfied in present experiments. Thus, a somewhat arbitrary 
separation of the observed scattering into l part due to critical fluctue :ions 
and a part due to short range order in the fluid is now required prior to a 
test of the predictions of universality (as defined above) . The best 
experiments to date seem to indicate significant differences between the 
predictions of latcice models and the behavior of real fluids (Warkulwiz et al., 
1974; Lin and Schmidt, 1974). Experiments with scattered light are now unable 
to satisfy the criterion mentioned above that the wavelength be much smaller 
than the correlation length. The correlation length increases rapidly as the 
critical point is reached. In a low gravity environment, the critical point 
could be approached much closer than is now possible; thus a much tighter test 
of this aspect of universality will be possible. 

Other questions which may be answered by experiments defined within this 
study pertain more specifically to aspects of fluid dynamics and structure. 

Hence the answers are less likely to have as broad an impact on condensed 
matter science than would a clear answer to the question of universality. 
Nevertheless, the properties of fluids are sufficiently alike at critical 
point > that questions about large classes of fluids may be answered with an 
experiment on one or two fluids. 

Questions of considerable interest are: What is the nature of the small 

anomalies which occur in the viscosity (Sengers, 1973), dielectric constant 



(Hartley et al,, 1974), refractive index (Hocken and Stell, 1973; Stell and 
lloye, 1974; Bedeaux and Mazur, 1973), and diameter of the coexistence curve 
(Weiner at n- . , 1974) (i.e. the average of the liquid and vapor densities) as 
the critical point is approached? There are theoretical and experimental 
controversies cn the nature of the anomaly for each of these properties, it 
is reasonable to believe that the more, definitive experimental results 
available from low-g experiments would have a significant impact on the 
theories for these properties in fluids. The experimental nature of the 
stronger critical poii. t anomalies in sound attenuation (Thoen and Garland, 

1974) and thermal conductivity (Sengers, 1973) is somewhat better understood 
from earth based experiments, hence, low-g environments are less likely to 
have an impact on the theory ,>f these properties. 

The dynamics of the process of macroscopic phase change in pure fluids is 
poorly understood in the vicinity 01 the critical pt it. A variety of 
questions remain to be answered. We wir 1 briefly mention a few. There is 
least one observation which suggests standard nucleation theories are 
fundamentally wrong near the critical point in pure fluids (Huang ejt al . , 1973). 
I 9 this pbservation correct? Other earthbound experiments are necessary in 
this area, but only in a low gravity environment will it be possible to obtain 
a homogeneous, macroscopic volume of a fluid under conditions sufficiently 
close to critical to yield a definitive study of this problem. The spinodel 
decomposition mechanism of phase separation occurs in alloys, glasses, and 
binary liquid mixtures (Schwartz et al., 1975). Does this process occur in 
pure fluids? If it does, it is ouite likely that tba study of spinodai 
decomposition in pure fluids (where many of the macroscopic parameters of the 
theory are well understood) will clarlf} the understanding of splnodal 
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decomposition in these other kinds of systems. Are there critical point 
anomalies in the process of bubble and droplet growth (or evaporation and 
condensation)? Do they influence the time of equilibration of macroscopic 
two-phase samples of pure fluids as has been suggested in two public . ions? 
(Dahl and Moldover, 1972; Brown and Meyei , 1972). What will be the dominant 
mechanisms for macroscopic phase sepa** ition upon cooling a pure fluid below 
the critical point if bouyant forces on bubbles and droplets are greatly 
reduced? Many other questions are possible in this poorly understood area 
of phaae separation. 

In fluid mixtures, important scientific questions analogous to the ones 
above exist. Specifically, the question: Do real fluid mixtures have the 

same "universal" thermodynamic and correlation function behavior as lattice 
models and/or pure fluids, is unanswered. This question is more difficult to 
answer for mixture*. h..:ause they have additional thermodynamic variables. For 
example the "quality" which disappears at the critical point may be thought of 
as being either a mass density difference or a composition difference between 
coexisting phases. In principle, if the critical point is approached closely 
eno ’h, either of these variables (or certain others) could be used to answer 
the question of universality. Questions analogous to those above may be asked 
about weak anomalies and transport properties in fluid mixtures. The role of 
gravity as an experimental limitation to answering these questions is much less 
clear in the case of mixtures that in the case of pure fluids. Tins will be 
discussed briefly below. 
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2. Opportunities Provided by a low-g Environment . 

2.1 Introduction 

In this section we will discuss in detail the manner In which gravity 
affects phase transition and critical phenomena experiments in fluids. We will 
see that gravity limits the closeness with which the critical point may be 
approached in all earthbound experiments. Thus a low-g environment will 
provide an opportunity for conducting experiments closer to the critical point 
than is possible on earth. 

Measurements -Z the equilibrium properties of pure fluids near the critical 
point can in principle encounter two distinct kinds of limitations because of 
the earth's gravitational field. One kind of limitation is essentially a 
technical one. All experiments measure average properties ove~ some finite 
height. Since the variation of fluid properties with height becomes 
increasingly large as the critical point is approached, this averaging causes 
increasing _.rors as the critical point is approached. In practical cases an 
important averaging error occurs even for optical measurements (which at first 
thought might be expected to average over a height of a wavelength of light 
~0.5 pm). The size of this "averaging error" depends both upon the property 
measured and the technique used to measure it. We will consider representative 
cases below. We will find that for nearly all experiments the closeness of 
approach to the critical point is limited by averaging errors. These errors 
may prevent the answering of many of the questions we raised in the last 
sec cion. 

A second kind of limitation to accuracy in the measurement of equilibrium 
properties of pure fluids which is imposed by the presence of the earth's 
gravitational field is an unavoidable or "intrinsic" limitation. Relatively 



simple considerations show that close enough to the critical point the 
correlation length (or the size of the density fluctuations) for the fluid is 
limited In the earth's gravity field much as it would be limited if the fluid 
were placed In a small container. This means that the properties of the fluid 
Itself are altered by the gravitational field. No Improvements in experimental 
technique will enable earthbound experiments to overcome this limitation. To 
date there are no experiments which are limited by this phenomenon, but it 
appears that this limit will be approached in the near future. 

Gravity enters into the study of the dynamics of phase changes in pure 
fluids and most fluid mixtures by causing relative motion of the two phases 
which are almost invariably of different densities. Thus in a 1-g environment 
a bubble or droplet will rapidly travel to the top or bottom of a macroscopic 
sample after It has grown to a size of the order of 10 pm. There are a 
variety of "levitation" techniques for studying droplet and bubble growth 
which circumvent sedimentation. Unfortunately, none have been demonstrated as 
being appropriate for studies of pure fluids near critical points where precise 
control of temperature and pressure are also required. 

If the study of binary liquid mixtures were conducted at equilibrium near 
their critical (consolute) points, "averaging errors" similar to those which 
occur in pure fluids, would appear. In practice, most binary mixture studies 
have been conducted at constant temperature but not in diffusive equilibrium 
so that a different kind of error is present. The diffusion constant tends 
towards zero at the critical point and diffusive equilibrium may take days or 
weeks in samples a 1 cm high moderately close to the critical point. It has 
recently been shown that the "pressure diffusion coefficient" diverges at the 
consolute point of binary liquid mixtures (Greer et al . , 1975; Giglio and 
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Vend r ami ni, 1975). This dlvsrgenc* Indicates "svsrsglog" srrors of « different 
kind may be encountered In liquid mixtures. 

Ve will now discuss these experimental limitations more quantitatively 
and in more detail. 


14 



2.2 Averaging Errors In "Bulk" Experiments . 

As a first and most simple illustration of averaging errors we now 
consider a hypothetical experiment to measure the density of a pure fluid as 
a function of height exactly at the critical temperature. Such a measurement 
together with the thermodynamic relation for the gravitation contribution to 
the chemical potential. 


du - - mg dz (2.2.1) 

(where g is the acceleration due to gravity, m is the molecular mass, z 
is the height coordinate Increasing upward, and u is the chemical potential 
per particle) would make it possible to test the predictions of many model 
equations of state. 

These model equations of state all predict (see appendix A) a critical 
ii 'therm of the form 


u-y. 


P D 
c 


lo-p. 


Sign (P-P c ) 


( 2 . 2 . 2 ) 


(Here u, P and p are the values of the chemical potential, pressure, and number 
density respectively, and the subscript c indicates the critical value; D is 
a numerical factor of order unity which varies from fluid to fluid and 5 is a 
"critical exponent" which will be the same for all fluids if universality holds for 
fluids in the critical region.) Fig. 2.1 shows the outcome of an ideal experiment o 
this kind with xenon. This figure was constructed using the numerical values for 
the parameters listed in appendix D. Note that the density can change by more 
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Fig. 2.1. The variation of the reduced density of xenon with height at the 
critical temperature. The calculated curve is displayed on two different scales. 
Note that near the critical density, a pressure change of 10~6 P will cause the 
density to change 5X and a pressure change of 10~9 P will cause C the density to 
change IX. c 
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than 5% as a result of a pressure change, AP, which is only one one-millionth 
of the critical pressure. At equilibrium this pressure change occurs (at 
normal earth gravity, g Q ) with a height change from 0.25 mm below the meniscus 
to 0,25 mm above the meniscus. We have chosen xenon for this numerical 
illc it ration (and others in this report) because it is one of the few fluids 
for which relatively complete data exist in the critical region. It is inert, 
available in high purity and has a critical temperature near room temperature. 
(COj and SF^ are other well studied fluids which are likely candidates for 
low-g experiments.) 

One practical instrument for the measurement of fluid densities in thermal 
equilibrium at g Q is a float densimeter (Greer et _al. , 1974). In this 
instrument the bouyant force on a float of known density is balanced with a 
spring or magnetic restoring force. The deflection of the spring or the 
current in the electromagnet is then a measure of p, the difference between the 
average float density and the average of the fluid density over the float's 
height. Quantitatively, 


/p fluid (z) o(z) dZ 
/o(z) dz 


(2.2.3) 


where o(z) is the cross sectional area of the float as a function of height. 
The outcome of a hypothetical experiment with a 1 mm high cylindrical float is 
ahown in Fig. 2.2 (left). The average density p (points in Fig. 2) differs 
substantially from the local density at the average height, 5, of the float as 
soon as any portion of the float overlaps the meniscus. An alternative 
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presentation of the "data" is of interest. We may think of this experiment as 
one which measures an average "susceptibility", x^» of the fluid. This 
susceptibility is thermodynamically related to the isothermal compressibility 
and plays a role in the study of fluids which is analogous to the role of the 
magnetic susceptibility near a magnetic critical point. By definition, 


* T 5 p2 *T " p2 (J T ) ( ^ T “ 0 T (2.2.4) 

In Fig. 2.2. (right) we compare XjU) with the average value of x computed from 
our hypothetical experimental data by 


X 


exp 



(2.2.5) 


It is Interesting to note that the experiment which measures density by 
averaging over a finite height is analogous to a spectroscopy experiment with 
a finite resolution. According to Fig. 2.2, the density experiment has only 
a limited ability to "resolve" the diverging susceptibility. Details smaller 
than the resolution are blurred. Yet, some information may be recovered on a 
scale smaller than the resolution If theoretical guidance is available. For 
example, if the shape and relative spacings of the spectral lines present in a 
spectral peak are known accurately in a*' ai\ce, their absolute spacing may be 
determined to much greater accuracy than they can be resolved by the instrument. 
Similarly, if Eq. (2.2.2) were known to describe the density profile and 
if techniques were of sufficient quality, accurate measurements of 6, D, and 
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P c could be obtained with a float of, say, 1 mm height. Unfortunately, the 
right hand side of Eq, (2,2,2) is only the first term in an expansion about 
the critical density. The higher order corrections to Eq, (2.2.2) are not 
known with certainty; thus, the interpretation of data of finite height 
resolution becomes subject to the uncertainties arising from different choices 
of correction terms. It is at this point that gravitational averaging makes 
experimental tests of universality difficult, 

A wide variety of realistic experiments have encountered the same density 
resolution problems we have illustrated here with the hypothetical density 
measuring expetlnent described above. For example Hohenberg and Barmatz [1972] 
have analyzed in detail the effects of gravity averaging upon measurements 
of the constant volume specific heat and of the low frequency velocity of sound. 
It turns out that both experiments are strongly affected by gravitational 
averaging at the critical temperature when the meniscus i9 within the 
experimental cell (see Fig. 2.3). The shortest suitable calorimeters 
constructed to date are about 1 mm high and the shortest low frequency 
velocity of sound resonator is under 4 ton high, so that both of these 
experiments will not resolve the density dependence of the quantity measured 
within about 4% of the critical density at normal earth gravity. 

The density resolution limit scales with sample height, h, and gravity 
according to the relation 



mg h p 1/6 




( 2 . 2 . 6 ) 


The scale factor for the density profile at the critical temperature, D P /(p mg), 

c c o 
o 

varies among 17 fluids (Appendix D) from 910 m in He to 10,400 m in H 2< For 
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In order to obtain C in J/moie K multiply tne oramate ny z.*. mis 
adapted from P.C. HoKenberg and M, Barmatz, Phys, Rev. A 6 , 289 (1972) 


xenon, it is 1310 m, a fairly low value among those fluids which are likely 
candidates for low g experiments. Because the exponent <5 is between 4 and 5 
it is clear that a substantial reduction in gh is required to improve density 
resolution at the critical temperature. Certain optical experiments 

(whose limitations we will discuss below) have averaging heights of micrometers 
rather than millimeters, thus optical experiments on earth become potentially 
competitive with bulk experiments when the latter are carried out in an 

-3 

environment of 10 g. 

We have illustrated the density resolution limits at the critical tempera- 
tures arising from vertical averaging. The actual resolution is limited at 
other temperatures as well. The shape of the resolution limited region is 
shown in a qualitative fashion in Fig. 2.4. The true form of the region, of 
course, depends upon the particular property being measured and upon the 
techniques used. It is straightforward to estimate the extent of the gravity 
affected region along the temperature axis as we have done for the density 
axis. We will do so in a qualitative' fashion here. Then we will indicate 
how a more precise criterion for averaging errors could be used to define the 
gravity affected region precisely. Such a precise criterion is formulated in 
Appendix C.3* 

A qualitative idea of the gravity excluded region can be obtained by 
noting that the scaling equations of state (Appendix A) indicate that 
asymptotically close to the critical point all thermodynamic quantities 
(which usually are functions of two variables, say temperature and density) 
may be expressed in the simple form r^f(0). Here r and 0 are parametric variables 
which are related to the temperature and density by nonlinear transformations. 
Roughly speaking, r is a measure of the distance from the critical point and 0 is 
a measure of the distance "around" the critical point from the coexistence curve. 


The 
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Fig. 2.4. Qualitative location of thermodynamic states In the temperature-density 
plane which are Inaccessible to bulk experiments because of gravitational 
averaging. Averaging errors for a 0.3 cm high sample are substantial beneath the 
curve marked r ■ 10"\ Averaging errors for a 0.03 cm high sample are substantial 
beneath the curve r - 10 “ 5. 



exponent p depends upon the property under consideration, and f(0) is an 
analytic function of 0 which also depends upon the property under consideration. 
In those cases (including the important quantities C v * and Xj) where p < 0 
and where f(0) is not a very strong function of 0, the r dependence of the 
properties becomes dominant near T c « It follows that divergent quantities will 
assume approximately the same value along a curve of constant r and such 
curves have been sketched in Fig. 2,4, It is also true that the gravity 
excluded region for many experiments is also approximated ! urve of 
constant r. We have already estimated the gravity limit fo. jrements of 

X T and C v in a sample at the critical temperature, thus locating point B 
on Fig. 2.4. To find the maximum temperature on the saim contour of 
constant r (point A on Fig. 2.4) we note that in one m del equation of state 
a contour of constant r is identical with a contour of constant x^ (Appendix 
A, Eq. A. 13). Thus in this model x T has identical values at B and A. (In 
other realistic models the value of x T will be nearly the same at B as at A) . 
Now the power law expression for x T on the line Ap* ■ 0 la 
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T 


£ 

P 



-Y 


(2.2.7) 


We equate this expression for x T with the expression for x T along the 
critical isotherm obtained by differentiating Eq. (2.2.2) 
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( 2 . 2 . 8 ) 



Using the values of the constants from the appendices we find that the temperature 

limit for bulk susceptibility measurements is roughly 0.028 K for a 1 mm high 

sample of xenon and that this limit scales with gravity and sample height at* 

(gh/g H )^" • If a contour of constant r is used as a gravity limit criterion 

O o 

other experiments will encounter gravity limits at some value of (T-T )/T which 

c c 

also scales as (gh/g H 

From the work of Hohenberg and Barmatz [1972] f the constant volume specific 

heat 01 a 1 cm high cylindrical sample of xenon begins to deviate from the 

specific heat of a zero height sample by several percent within 50 millikelvin 

of the critical temperature; for a 1 mm high sample significant error will 

occur within 12 mK of T c# Naturally the exact temperature depends on the 

criterion of accuracy, but with a fixed criterion of accuracy the temperature 

0 

of closest permissible approach to T will scale as (gh) ‘ Thus a 

c 

substantial improvement in temperature resolution possible bg doing C 
experiments in a lcw-g environment . In particular, the determination of in 
specific heat experiments would be greatly facilitated in a low g environment. 
Instead of approximating the region of severe gravitational averaging by 
a curve of constant r (es we have just done), a precise calculation could have 
been done. The calculation would involve the following steps: 1) Identify 

the quantity to be measured (say' Q) and the precision with whi~h the quantity 
is to be known; frequently the desired precision may be expressed as a small 
fraction, p, of the quantity Q, 2) Examine the measurement technique to 
determine over what range of heights it averages (say 0 £ z £ h) and with what 
function (say w(z)) it weights measurements at each height. 3) Compute the 
values of p and T for which the inequality 


1 

h 



h 

w(z) Q (p(z), T) dz - Q (p(h/2) , T) £ p Q(p(h/2), T) 


o 


(2.2.9) 
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is satisfied. Here we have assuned that the weighted average of Q which the 
experiment measures will be assigned to the thermodynamic state (p(h/2), T) 
at the midheight of the experiment. Exactly this calculation is carried out 
in Appendix C where Q is taken to be the correlation length and w(z) is 
assumed to be unity. A somewhat more complex calculation id required if Q is 
a quantity measured by experimentally taking a temperature derivative (such 
as the constant volume specific heat). This problem is discussed by Hohenberg 

Barmatz [1972]. 

We have Just completed several illustrations of how gravitational 
averagin} 'ufluences the measurement of thermodynamic quantities. It is 
important to note that most measurements of transport properties will also 
be subject to limited density and temperature resolutions in the 
earth's gravity. For example the density dependence of the viscosity or of 
the thermal conductivity cannot be measured within 4-5% of P c> if the viscosity 
or conductivity apparatus is 1 mm high at earth-normal gravity. It is also true 
that the density dependence of the turbidity (or total light scattering 
.ntenslty) cannot be measured within a few percent of p £ if the scattering 
volume is 1 mm high (Leung and Miller, 1975) because of the vertical averaging 
that occurs in the scattering volume. On the other hand quite different 
limitations apply to the measurement of thermodynamic and transport properties 
by optical techniques with very fine spatial resolution. They are 
discussed in the subsequent section. 





2,3 Limitations on Optical Experiments due to Gravity Induced Refractive 

Index Gradients. 

As a fluid is brought near its critical points, its coefficient of 
Isothermal compressibility diverges and, in a gravitational field, the fluid 
becomes compressed under its own weight. The density gradients thus produced 
place strong limitations on the validity of measurements, even with optical 
probes, very near the critical point. The most serious limitation may be 
termed a "thick cell effect" and applies to all optical experiments including 
measurements of the phase, intensity and spectral characteristics of light 
passed through the fluid. Light which is directed horizontally into an 
optical cell filled with fluid near the critical point will be deflected down- 
ward by the index of refraction gradient (resulting from the density gradient). 
The angle of deflection is proportional to the thickness of the fluid layer in 
the cell. This effect is precisely the same one which enables the sun to be 
seen above the horizon several minutes after it has "set" according to 
astronomical calculations. As the critical point is approached, the density 
gradients and the deflections become so large that light "rays" pass through 
layers of fluid of widely varying density. Then it is no longer possible to 
relate the Intensity, the spectrum, or the phase of the light emergent from the 
cell to the density or any other local thermodynamic variable. The measurable 
optical properties of the fluid become complex gravitational averages. In this 
section we discuss and quantify these limitations for optical experiments 
using a combination of analytical and numerical techniques. 

Consider a sample of a dielectric fluid, in a cell with plane parallel 
optical windows, with an average density, p q , closely approximating the 
critical density. When the sample temperature is held constant near T^, the 
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fluid is compressed under its own weight and an equilibrium density distribution, 
p(z,T), is formed. Since, in a gravitational field, the reduced chemical 
potential u is directly proportional to z, the density profile is directly 
given by the equation of state p(y,T). Furthermore, if current theory is 
correct (Appendix A) such profiles are antisymmetric about some density 

p (p is not necessarily the critical density and in general will be a 

cen cert 

function of the temperature. For this analysis we will assume that p * p ) 

cen c 

which occurs at the height where the compressibility has a maximum. The plane 
at this height, which we call the centrus, will serve as the origin for our 
vertical coordinate system. In Figure 2.5 we represent such a profile 
schematically for some supercritical temperature T > T^ (All of our analysis 
will deal only with the one phase region.) Because the refractive index of the 
fluid is related to its density by the Lorentz-Lorenz (Lorenz, 1880 ; Lorentz, 
1952; Hocken and Stell, 1973 ; Larsen et al., 1965 ) formula, there is a 
corresponding refractive index profile n(z,t) in the sample. Now we 
illuminate the optical cell with a monochromatic horizontal (normal to the z 
axis) plane wave of light, and examine the trajectories of rays through the 
profile at various heights. 

In Fig. 2.6 we illustrate this situation with a simple ray picture. A 
ray enters the cell horizontally at height z’ above the centrus and is 
refracted downward on a curving trajectory until it emerges at a height z. 

The output angle of the real ray will be proportional to the average refractive 
index gradient ’seen 1 by the ray and Its phase shift will be determined by an 
average over its optical path. Also shown in the same figure are two ideal 
rays which we call ’thin cell’ rays. Both are equivalent to the real ray when 
the cell is infinitely thin or the gradient so small that z’ z. They are 
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Fig. 2.6. Rays through a density profile. "Real ray" denotes the path 
traversed by a light ray incident horizontally upon an optical cell filled 
with fluid near the critical point. The other rays are defined as follows: the 
angle equivalent ray (TER) is that thin cell ray which is refracted at an angle 
equal to the real ray. It traverses the cell horizontally at a height (z + A z) 
where the local refracti\e index gradient is equal to the average gradient seen 
by the real ray; the height equivalent (ZER) is that thin cell ray which traverses 
the cell horizontally at height z. 
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defined as follows: the "angle equivalent*’ ray (TER) is the thin cell ray which 

is refracted at an angle equal to the real ray; it traverses the cell 
horizontally at a height (z + Az) where the local refractive index gradient is 
equal to the average gradient seen by the real ray . The '’height equivalent*' ray 
(ZER) is that thin cell ray which traverses the cell horizontally at height z. 

These thin cell rays are introduced for courarison purposes since for them 
the relation of optical observables to thermodynamic quantities is particularly 
simple (Estler et _al. , 1975). For instance, if we expand the Lorentz-Lorenz 
formula, we find to a good approximation that 


n 


n c * n l 



* n i Ap* + . . . 


(n 2 - 1) (n 2 + 2) 

(n a — ^ ). 

6 n 

c 


Then the index gradient becomes 


( 2 . 3 . 1 ) 


n 

dn do * 1 

dz * “1 dT~ * IT 

c 




(2.3.2) 


where u* is the reduced chemical potential, Xj* the reduced susceptibility 

and H q * P c /(p c rag), the scale height for the fluid being studied. In the thin cell 

limit the ray ZER simply comes out at height z with an angle proportional to 

A 

X T at that height, i.e. 
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X *(z) • for the ZER, (2.3.3) 

T dn^ 

here d is the cell thickness. 

Fig. 2.7 offers a comparison between angles of rays as a function of z 

determined from equation (2.3.3) and angles determined numerically by 

tracing rays (Born and Wolf, 1975) through the same density profiles. The 

profiles were generated using the Stony Brook parametric equation of state (Wilcox 

and Estler, 1971) with parameters obtained from actual profile measurements on xenon 

(Estler et al. , 1975). The parameters used are given in Appendix E. The 

T-T _ 6 

computation was performed with AT* ■ — — - m 5x10 . The solid line represents 

c 

the thin cell results for the ZER (Eq. 2.3.3) and terminates in a rounded point 
an angle of about 17 radians, that is, off the graph. The dashed line is 
the result of the ray tracing computation. The crosshatched region contains 
rays that enter the cell close to the centrus above and below and emerge at 
approximately the same height but at very different angles; the observer sees 
a bright band at this height with a darker band abo^e it. In this region the 
gradient is so large that the rays are simply bent right out of the dark region 
and into the bright band. 

Fig. 2.8 shows the same results for the optical phase as a function of 
height. Again the solid line is the thin cell result, i.e. neglecting any beam 
bending effects, (cf. (2.3.3) and the dashed line is the result of ray tracing. 

The optical phase relative to that at the centrus (A$) is simply related to 
the density in the thin cell approximation: 

M 1 * k d Oj p*(z) (2.3.4) 
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Fig. 2.7. The angles of emergent rays as a function of height. The solid 
line represents the thin cell results for the ZER (eq. 2.3.3) and terminates in 
a rounded point at an angle of about .17 radians, that is, off the graph. The 
dashed line is the result of the ray tracing computation. The crosshatched 
region contains rays that enter the cell close to the centrus above and below 
and emerge at approximately the same height hut at very different angles. 




Fig. 2.8. Optical phase as a fur tlon of height. The phase profile calculated 
In the thin cell approximation (solid curve) has the same shape as the density 
vs height profile. The observed profile Is represented by the dashed curve. 

The crosshatched area corresponds to the crosshatched area In Fig. 2.7. 
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Hence, the solid curve Is, except for sign and units, simply the existing 
density profile. The dashed line is then the apparent density profile deduced 
from actual measurements and in fact such "kinked" profiles have been reported 
in the literature (Lorentzen and .Hansen, 1966). 

Figures 2.7 and 2.8 offer a good visual picture of what happens to light 
traversing a fluid near its critical point. They do not, however, offer 
simple quantitative guidelines as to which regions of the thermodynamic 
space are inaccessible to the earth bound experimenter. To provide these 
guidelines we found it convenient to compare the real ray with the angle 
equivalent ray (TER) rather than to the ZER as in Figs. 2.7 and 2.8. 

As our computations progressed we discovered numerically a simple 
relationship between Az (the height difference between the TER and the 
real ray, see Fig. 2.6) and the output angle of the real ray (or TER since 
they are defined to be the same)* We then derived thlB expression for Az by 
solving the ray tracing equations analytically in a medium where dn/dz was a 
slowly varying function of height. We found that if 



then 


Az iv 


,2 

V * 
3HnJ X r 


(2.3.6) 


As long as equation (2.3.6) holds, the real ray carries average fluid 
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information equal to the local information at a height Az above the height of 
the real ray's emergence. This equation la valid both above and below the 
centrus. 

We then performed computer experiments to test the domain of validity of 
(2.3.6) using numerical parameters appropriate for xenon. We found that (2.3.6) 
breaks down when Az reaches a certain value nearly Independent of the temperature 
and the fluid studied. The results are plotted in Fig. 2.9 for xenon. In this 
figure the solid dark line represents the Az that would be calculated from 
(2.3.6) using the fluid's maximum compressibility. That is 


Az 


max 


V . 

3 “"c 


n. d -y 

Tita • 

c 


(2.3.7) 


The region above and to the right of this curve is nonphysical: no real 
thermodynamic states exist in this region. The area below this curve and to the 
left is physical but mostly inaccessible to the earth bound experimenter. On 
this plot the critical point is at infinite Az. The dashed line at Az & .002 cm 
represents the contour where (2.3.6) becomes incorrect by IX. Above this line 
the errors grow very rapidly. The region below this line is experimentally 
accessible and hers the properties of the interference pattern may be simply 
related to local thermodynamic properties of the fluid. As a simple rule of thumb 
we find that the observed interference pattern can no longer be simply related 
to local thermodynamic fluid properties when the reduced compressibility 
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Fig. 2.9. Thermodynamic states acce. c 4 .bl» to optical experiments. The 
states with sufficiently small values ot the scaled compressibility x^ are 
accessible. Large values of x<j. <1° not occur to the upper right of the solid 
curve (i.e. far from the critical point). 
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becomes larger than 


* r 

X - (.002 cm) 1 . (2,3.1) 

max n^d 

Ic is convenient to tabulate some examples for comparison with gravitational 


limits provided elsewhere in this 

* 

report We do 

so here for 
* 

xenon 
4 * 

Cel- Thickness 

x T 

0 max 

^ p min 

AT . 

min 

cm. 

max 

milliradlans 

(At* - 0) 

(Ap* - 0) 

1.0 

2.6xl0 3 

22 

0.0471 

1.6xl0 -4 

0.1 

2.6xl0 5 

220 

0.0119 

3.3xl0~ 6 


Thus any optic*., experiment using a 1 cm thick cell is subject to limitations 
comparable to a bulk experiment in a 0.1 cm high cell at g Q . An optical 
experiment in a 0 1 cm thick cell is subject to limitations comparable to a 

_ T 

bulk experiment in a 3 ym high cell at g Q or e 3 mm high cell at 10 g Q » 

Optical cells with properties suitable for critical point experiments (Hocken 
et al. , 1975) have been fabricated with a thickness of 3 mm. It seems likely 
that 1 ran thick cells could be made. The feasibility of making much thinner 
cells should be examined. One may expect that fluid samples in extremely thin 
cells would be subject to contamination from the cell walls because of the large 
surface to volume ratio. 
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2,U Light Scattering Measurements and Some Gravity Related Limitations . 

Scattering measurements provide a means to measure the space-time 
dependence of fluctuations. Since much of the fundamental picture of critical 
fluctuation phenomena deals with the distance scale over which fluctuations 
are correlated (space) and the tine scale over which fluctuations grow and 
decay (time) *scattering spectroscopy continue! to be an important tool for 
studying critical phenomena. 

Near the gas-liquid critical point a fluid exhibits large fluctuations in 
the density. The magnitude of these fluctuations is proportional to the 
compressibility. The spatial extent of these fluctuations can be characterized 
by a correlation angth £, As shown in Appendix C this correlf.ti.on length is 
related to the compressibility by 


c - £ (r _1 x!) 2 ' n (2.4.1) 

O I 

* 

where is the dimensionless symmetrized compressibility defined in A. 2. The 

constants and T are the amplitudes in the power laws (C.4) and (A. 4c) for the 

correlation length and compressibility along the critical isochore. The constant 

£ is a microscopic distance of the order of the range of the intermolecular forces, 
o 

The spectral intensity as a function of scattering angle and frequency is 
proportional to the spatial and temporal Fourier transform of the time 
dependent correlation function of the order parameter. In a one component 
fluid the integrated intensity of the scattering is proportional to the 


static structure factor 


(2.4.2) 


i 


S(k) - I dR e ikR G(R) 

where the correlation function G(R) Is defined In Appendix C. The wave 
number k Is related to the scattering angle 6 by 

k - 2 k Q sin | (2.4.3) 

where k is the wave number of the incident radiation, 
o 

In principle S(k) is a function of temperature AT*, density Ap* and wave 
number k. However, according to the scaling hypothesis the structure factor 
near the critical point can be written in the form (Fisher, 1967) 

S(k) - C 2 ' n <KkO (2.4.4) 

where £ is the correlation length. This reduction of the description is vaxid 

for any thermodynamic pat;, in the AT* - Ap* plane of constant scaling variable, 

1/0 

i.e. AT* proportional to | Ap * | . Here we shall consider specifically the 

properties of scattered light at ~he critical isochore Ap* - 0. 

The scaled c unction <Kk£> is known theoretically (Tracy and McCoy, 19 7 5) 
for the 2 dimensional Ising model for all values of k£. However, the 
mathematical form of this function for fluids is not known except for the fact 
tnat :s must approach a constant in the limit k£ -► 0 and that it must vary as 
(kO~ 2+ ’’ for kC ->• ». 

It has turned out to be very difficult to determine a definitive value for 
the exponet t n. A precise knowledge of this exponent has a direct bearing 
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on the question of universality of the pair-correlation function. 

In order to datsmlne this experimental n with any accuracy, measurements 
are needed for kC sufficiently large eo that the function <Kkt) can be 
approximated by its asymptotic behavior (Tracy and McCoy, 1975). Due to 
limitations caused by multiple scattering and the gravity Induced density 
gradients, it has been impossible to satisfy thxS condition for earch bound 
light scattering experiments. It has been possible to reach sufficiently 
large values of k with neutron scattering (Warkulvlz et al. , 1975); however, 
here one has difficulties in satisfying the condition that k * must be large 
compared to the range of the intermolecular forces. 

Of the several available scattering techniques (light. X-rays and neutron 
scattering), optical light scattering using laser sources seems to be most 
feasible one for space experiments. It is the only scattering technique with 
sufficient spectral resolution to examine the narrow spectrum of the 
fluctuations, can be done with easily portable sources (lasers), requires no 
radiation shielding, and uses windows compatible with other optical records 
desired of the sample behavior. The principle drawback to light scattering 
mef- .rements is multiple scattering: the extinct' in of the beam gives a severe 
lower limit on the AT* which can be reached. Numerical estimates for this 
limit are presented later in this section. 

Using optical beating techniques, one can measure the spectrum of the 
quasi-elastic Raylel^' line (Cummins and Swinney, 1970). The width of this 
line is related to the diffusive decay constant of the fluctuation of the order 
parsu-iter. A survey of the current status of the experimental work on 
Rayleigh scattering near the critical point of fluids has been given by 
Swinney and Henry [1973], 
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In interpretating such experiments one usually introduces the following 
assumptions. First it is assumed that fluctuations still decay 
exponentially in time so that the dynamic structure factor has ..te form 


2 r 

S(k»G>) = S(k) - 5 — V" (2.4.5) 

r + w 

S 

where T * is the decay rate of the entropy fluctuations. (Svinney and Henry, 1973). 
Secondly, one uses the assumption of dynamical scaling to write the decay rate 

T in the scaled form 

s 

r = r 3 * r (kO ’ ( 2 .A. 6 ) 

It should be pointed out that these assumptions are not valid rigorously, and 
that one expects to see deviations if the critical point is approached 
sufficiently closely. Theories for the scaled function ♦p(kS) have been 
developed by Kawasaki and coworkers and by Ferrel and Perl; for a survey of the 
literature the reader is referred to the article by Swinney and Henry [1973] « 

The various theories differ in their conclusions, but it has been impossible 
to discriminate between the theories on the basis of earthbound experiments. 

For light scattering measurements we need to consider the limitations on 
the attainable experimental precision due to gradients in the density Ap*, 
turbidity of the sample, and gradients in the correlation length £. Our 
purpose here is to estimate how much the accessible range in AT* and Ap* can 
be extended by conducting light scattering experiments at reduced g. 
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We first consider the variation in density. In a light scattering 
experiment close to T c one can consider using a weakly focused beam with a 
diameter of 2 100 urn. Stronger focusing would require much lower power levels 
to avoid local heating particularly at window surfaces. Thus there exists a 
practical height or spatial scale limit for critical point light scattering 
measurements of about. 100 ps. 

The resolution of light scattering experiments has been exploited in earth- 
bound experiments by moving the beam as a function of height, thus obtaining 
local equilibrium isothermal scattering Intensities (White and Maccabee, 1975) 
and spectra (Swinney and Henry, 1973; and Rim, et al. . 1974). In order for 
the density change over a height to be within a precision p of the average 
density in that region for a sample in a gravitation acceleration g*(g* * g/g Q ) , 
we require 


d Ap* 
dz 


h £ p 


(2.4.7) 


Using the restricted cubic model and Eqs. A. 13 and B.7 of the Appendices, we 
find 


ic - 


( P W 


(2.4.8) 
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For xenon , using parameters of Appendix D, ve find 


r * 6.2xl0' 4 (JlSl)- 84 (2.4.9) 

with h in meters. Thus this density limit is a contour of constant r In the 
(AT*, Ap*) plane as discussed in Section 2.2. The following 

table gives examples of the lower bounds to r and corresponding AT* * r (with 

-4 

Ap* * 0) and Ap* (with AT* * 0), for p * .01 = IX, h ■ 10 m = 100pm corres- 
ponding to precision scattering experiments on earth (lg), worst case Shuttle 

—3 —6 

accelerations (10 g) and best case Shuttle accelerations (10 g). 


Table 2.4.1 


g* " g/g„ 

r - AT* 

Ap* 

o 

(Ap* * 0) 

(AT* - 0) 

1 

1.3xlO~ 5 

1.9xl0~ 2 

10~ 3 

3.9xl0 -8 

2.4xl0 -3 

M 

O 

1 

ON 

1.2xl0~ 10 

3.0X10" 4 


Next we examine the effect of the strong (critical opalescense) scattering 
of light in the sample. The density fluctuation scattering cross section per 
unit volume per stearadlan of a fluid is given by 


nv») 

dr 


= *2 


$1n 2 4> /, 3e\2 


k B T K T 


(2.4.10) 
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where X is the wavelength of the incident beam (in vacuum), e is the optical 
dielectric constant at X and is the isothermal compressibility and ■$ the 
angle between the incident light polarization and the scattered wavevector. 

Near the critical point has a strong temperature dependence leading to 
the result that the scattering intensity diverges like K^. For k » 2tt/X << £ \ 
this dipole cross section can be readily integrated over all angles to give the 
scattering attenuation coefficient called the turbidity, t 

' ■ *r ^ (0 ii >: k s T *1 < 2 -*- 11 > 


For the present estimates we have neglected the angle dependence of the 
scattering due to the correlation length £ (see Tuglielli and Fore, 1970 and 
Cannell, 1975). 

Like the cross section, the turbidity follows the temperature dependence 
of K^. Using the cubic model (Appendix A. 3) we can write for K^. 






(2.4.12) 


Since close to the critical point p 3 P c> we find that t is proportional to r r . 

In light scattering studies the intensity corrections due to turbidity 
become severe for a path length l such that £ 1. At approximately the same 
t£ limit the corrections due to multiply scattered light reaching the detector 
become severe. Thus we see that scattering experiments are limited to studying 
the samples for AT* - Ap* such that t < i 

We have evaluated the turbidity and the resulting limitations tor the case of 
xenon using the parameters of Appendix E at A = 6328A (He-Ne laser), and 
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and approximating (p— ) (e -1). We find 

op C 

T - 2. 51x1 0~ 3 r" 1 - 19 (m* 1 ) . (2.4.13) 

For an optical path t (m) we conclude that r must satisfy 

BA 

r * 6.53x10 £* . (2.4.14) 

-2 

Typical optical cells at present have 1-10 m (1 cm). More specialized cells 

_3 

have been made with £ * 10 m (1 ran). It may be feasible to work with 
L - 10 (100pm). With such a thin cell one must carefully discriminate against 

scattering from the inner cell walls. These distances have been used to 
calculate the lower bounds of r and corresf r nding AT* and Ap*. For xenon and 
X - 6328a we find: 


Table 2.4.2 


Turbidity approach bounds at A » 

6328A in Xenon 

£(m) 

r * AT* 

Ap* 


(Ap - 0) 

(AT* - 0) 

io -2 

1.4 xl0~ 4 

4.4xl0~ 2 

io -3 

2.0xl0* 5 

2.2xl0 -2 

io ' 4 

2.8xl0 -6 

l.lxl0" 2) 
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Comparison of Tables 2.4.1 and 2.4.2 ehovs that at lg, aim optical 
path cell of xenon has both the density profile and turbidity limit of about 
r - l.xlO \ This limit corresponds to AT » 3 mK. In order to take advantage 

-3 

of decreased density gradients at 10 g one would need a cell with optical path 
-4 

less than 100pm = 10 m. The ability to fabricate such a cell has not been 
demonstrated. Thus we see that turbidity and multiple scattering place a severe 
bound to scattering experiments closer to the critical point. On earth at lg this 
has not been as severe because the sample density gradients have caused most of 
the sample to be off the critical density giving reduced scattering and 
turbidity. In a low g environment one expects to see the sample uniformly 
opalescent and becoming quite opaque to visible wavelengths as the critical 
point is approached. The fact that the sample is uniformly opalescent at low g 
should make the analysis of the influence of multiple scattering much more 
manageable. The new theoretical work on multiple scattering (Reith and Swinney, 

1975; Bray and Chang, 1975) will be helpful in the analysis of low g scattering 
data. 

It Is clear that for studies close to the critical point one needs to 

reduce the observed turbidity. This can be done with the choice of X or choice 

-4 

of fluid with a small defractive index. Since the turbidity varies as X , the 

turbidity at a given AT* could be greatly decreased by using longer wavelengths. 

Small continuous He - Ne lasers exist at 1.152pm and 3.391pm. Photomultipliers 

still work at. 1pm so this choice would be favored for spectral measurements and 
would decrease the tubidity by 10. 9x and the r bound by 7.5x. At 3pm photoconduct ive 
detectors would be required. The quantum efficiency may be so low that beating 
spectra of the critical fluctuations may not be possible. Also the choise of window 
materials becomes more restricted. The following table shows the turbidity 
bounds of AT*, Ap* in xenon: for electromagnetic radiation with a wavelength of 
3pm. 
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Table 2.4.3 


Turbidity approach bounds at X * 3.391v in xenon 


i(m) 

r - AT* 

Ap* 


(Ap* -0) 

(AT -0) 

lO" 2 

4. 8*10" 7 

6.0xl(f 3 

Ilf 3 

7. OxlO -8 

3.0xl0“ 3 

1 

O 

i-J 

l.OxlO" 8 

1.5xlO" 3 


-3 

Table 2.4.3 Indicates that at 3.391ua we could exploit the 10 g environment 

—5 —8 —2 

with lOOym path cells, over a new range of AT* • 10 to 10 and Ap* - 2x10 

• —6 
to 2x10 . With 6328A we could only go in to AT* - 10 . 

Alternatively we can consider going to a fluid wiu.i lower dielectric 
2 

constant, e - n . The best choice (lowest e) would be the rare isotope of 
helium, He 3 . Here n^ 2 g * w 1.0108 in the fluid near T c « (We have taken n £ for 
He 4 from the data of Edwards and Woodbury [1963] and scaled it for He J using the 
Clausius- Mosotti relation with the critical densities from Appendix D) . Thus at 
6328 a we expect T H€ / T Xe * (e^ -1) /(e^-l) “ One could further reduce by 

moving out to lpm or 3pm. Thus, the use of helium would help overcome the 
increase in turbidity encountered upon approaching the critical point at the 
cost of doing experiments at cryogenic temperatures (3.3 K). 

The angle dependence of scattering intensity and linewidth closer to the # 
critical point are of interest because of questions concerning the theory of 
the fluctuations in the non-hydrodynamic region. The angle dependence of 
intensity is written in scaled form in (2.4.4), At a fixed AT* f Ap* the angle 
dependence arises from the function <j(x) with x ■ k£ and k ■ 2 k Q sin 0/2. 

Similarly, the scattering angle dependence of the decay rate i9 given 
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by (2.4.6) and comes from x in the function <|> r (x). If the measurements are 
to be used to determine £ or to test the form of <J»(x) or i^(x) it is essential 

that we know x precisely. If £ varies over the sample or specifically over the 

diameter of the probing beam in the sample due to gravity gradients then x is 

correspondingly poorly defined. Assuming for the moment that k is well 

defined through suitable coll imat ion of the Incident and scattered beams, the 

measurement precision of x • k£ will be limited by the uniformity of £. 

-4 

We ask that £ be constant to a precision p over a height h ■ 10 m = 100um 
or 


h S P 4 


(2.4.15) 


Using the results oi Appendix C for C~^) T one can use (2. A. 15) to give 

r(0) and then transfer to AT* f Ap*, We have done this for xenon with p - 0.1, 
-A 

h * 10 m giving the lower bound AT*, Ap* contours shown in Fig, 2.10 for 
—3 —6 

g* ■ 1, 10 and 10 . Comparing the figure to Table 2.4.3 shows that even 

if the turbidity limit is lowered by changing to 3um wavelengths we will 

_3 

begin to be limited for g* ■ 10 by the correlation length gradient. 

There is a further effect of turbidity in determining x. This occurs 
because k is well defined only for non-attenuated Incident and scattered 
waves. Attenuated waves give a spread to k. To have k defined to 1% we 
need k > 10 t. Since k = x/£, we require 

* A * <^> c. (2.4.16) 
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Experiments with a lOOum diameter beam in the region below and to the right 
of each contour are severely affected by averaging over correlation lengths 
at the values of g* indicated. 
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This Is lesst restrictive for large x. The condition (2.4,16) gives a lower 
bound for r which is converted into AT* for Ap* • 0 and Ap* for AT* » 0. 

Table 2,4.4 gives the results for X » 6328 a and X - 3.391pm for the case x » 1. 
Note they are independent of g. 

Table 2.4.4 

Scattering vector approach bounds 


X 

r ■ AT* 

Ap* 


(Ap* » 0) 

(AT* • 0) 

6328A 

6. 3x10“ 7 

6. 6x10" 1 * 3 

3pm 

1.6x10“® 

1.8x10" 3 


The results in Table 2.4.4 indicate that the change in wavelength from 

6328 to 3pm would keep the smearing of k to an acceptable limit down to the 

-4 

turbidity /multiple scattering limit for t of the order or 10 ' m or 10 m. 

We can summarize the results reported here by stating that with the most favorable 

geometry of 100pm path cells and 3um radiat' n light scattering experiments could 

3 3 

approach closer to the critical point by 10 in r corresponding to 10 in AT* 

1 -3 

and 10 in Ap* in a 10 g environment. At present , experimental limitations 
would prevent us from fully utilizing the minimum accelerations of 10 8 g. Even 

at 3pm a requirements remains from j* i to achieve the minimum r allowed by 
density and turbidity we muet f ill the cell accurately to p c to within .02%. 

This is certainly possible but requites more care than does sample loading for 
earth bound experiments. 
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3. Limitations in critical-region experiments due to modifications of 

fluid properties by a gravitational field . 

In the preceding sections we have commented on the technical complications 
that are encountered in various experimental methods near the critical 
point of a fluid in a gravitational field. However, in addition, 
intrinsic limitations in earth-bound experiments exist due to the fact 
that the gravitational field modifies the fluid properties in the immediate 
vicinity of the critical point. 

As mentioned in Section 2.4, upon approaching the critical point the 
increase of the compressibility is accompanied by an increase in the 
size of the fluctuations that extend over a correlation length £. If 
the system were homogeneous and in true thermodynamic equilibrium the 
compressibility and, hence, the correlation length would actually diverge 
at the critical point. However, the presence of the gravitational field 
prevents the fluctuations from growing indefinitely and compressibility 
and correlation length will In fact remain finite. Thus the presence of 
the gravitational field causes round off effects which change the 
nature of the thermodynamic behavior in the immediate vicinity of the 
critical point. It is the purpose of this Section to estimate the range 
in temperature and density where experiments will be affected by these 
round off phenomena and to elucidate how this range depends on the 
magnitude of the gravitational field. 

As explained in Section 2.2 the gravitational field induces a 
density gradient in a fluid near the critical point. When the density 
does not change too rapidly as a function of height, one may assume that 
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the local thermodynamic properties of the fluid at a given Jevei are the 
same as that of a macroscopic homogeneous system with the sare values of 
density and temperature. Under these conditions the density profile is 
determined by the equations given in Appendix 5 and measurements of the 
local fluid properties in a gravitational field does provide information 
on the thermodynamic behavior of a homogeneous fluid in the absence of 
a gravitational field. 

The assumption of local thermodynamic equilibrium is justified when 

the fluid is locally homogeneous over distances of the order of the 

correlation length, but will break down when the fluid properties begin 

to vary over distances of the order of the correlation length. In 

Fig. 3.1 we show calculated density profiles for xenon in the earth 1 s 

gravitational field at three reduced temperatures, namexy aT* = 10 , 

_6 

10 and 10 . The existence of a density gradient implies that also the 

correlation length £ varies as a function of the height. In Fig. 3.2 we 
show the correlation length £ as a function of the height for xenon in 
the earth’s gravitational field. The curves are calculated in the 
approximation that the assumption of local thermodynamic equilibrium 
remains valid. However, in the dashed part of the curves the correlation 
length varies so rapidly that the system can no longer assume local 
equilibrium states that are homogeneous over the range of £. Under 
these conditions the laws of thermodynamics no longer suffice to specify 
the nature of the equilibrium states. The complications that arise when 
the macroscopic thermodynamic relations can no longer be applied at the 
local level, will be referred to as non local effects. 

Such effects will enter when the correlation length £ starts to 
vary over its own height h ■ £. It thus follows from (C.19) that for 
measurements conducted with precision p non local effects will be encountered 


when 
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Fig. 3.1 Calculated density profiles for xenon In the earth's gravitational 
field and assuming local thermodynamic equilibrium. In the dashed 
part of the curves the density varies so rapidly that the assumption 
of local thermodynamic equilibrium no longer applies. 







(3.1) 


X 

j 

t 

I 


I 



a P 


In terms of the parameters r and 0 of the restricted cubic model 
equation of state, it follows from (C.16) that non local effects can be 
avoided only if 


r v+0<S 



[ 2b 2 *e 2 

[266b 2 2 0 2 O-0 2 )+(l-30 2 )(l-b 2 2 0 2) 




(3. 


Using the xenon parameters given in Appendix E, we conclude that at the 
critical temperature, AT* - 0 (0 » *\> 2 *), non local effects will be 

encountered unless 


| Ap*| * 0.007 



_£ 

v+86 


0 007 (p“) 


0.163 


(3.3) 


In Fig. 3.3 we indicate the region in the AT* - Ap* plane where the 
behavior of the fluid is modified by non local effects and, hence, where 
its properties will be fundamentally different from a fluid in the 
absence of gravity. 

Near the critical isochore, Ap* » 0, condition (3.1) becomes 
unrealistic because of the rapid variation of 8£/3z and we should, 
instead, consider the Integrated form 



The distance l.z over which the change in £ is smaller than p is 
P 

derived in Appendix C. In particular we require 
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0 0.01 0.02 0.03 


( P • P c ) / Pc 

Fig. 3.3 Region in the temperature-density plane where the fluid properties 

are modified by non local effects for g* = 1 (earth's gravitational 

field) and g* = 10”^ The curves refer* to experiments for xenon 

_2 

with a precision at the 1% level (p = 10 ). 
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(3.6) 


' 'aT*I" w 2.\i 

•o' • p 


where a* is defined in (C.13). For xenon we obtain 
P 


I AT* | £ 0.4x1 0~ 6 (2V +SS 

•i> 


= 0.4x1 0" 6 (S^)°« 46 


(3.7) 


This equation determines the intercept of the contourt in Fig. 3.3 with 
the aT* axis. 

It follows from (C.lla) and (3.2) that the system will only satisfy 
the conditions for local thermodynamic equilibrium when 


4 s 1 . 5 xl 0 -6 (|*) v+66 m = 1 . 5 xl 0" 6 ( jj *)°‘ 29 * ( 3 * 8 ) 

When the critical point is approached more closely, the increase of the 
range of the fluctuations will be suppressed by non local effects. It 
seems plausible to assume that the correlation length £ in a fluid in 
a gravitational field cannot grow any further when (H/3z)^ - 1 in the 
local equilibrium approximation. Hence, in a fluid in a gravitational 
field, the correlation length will always remain finite and of the order 
of 


5 < 1 .5x1 O’ 6 g*' 0 ' 29 m 


(3.9) 
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It should be noted that the maximum correlation length attainable will 
only increase inversely proportional to the cube root of the gravitational 
field. 

In this section we have considered the intrinsic limitations of 
thermodynamic experiments near the critical point of a fluid in a 
gravitational field. At present * more stringent limitations are imposed 
by the complications associated with the state of the art of the various 
available experimental methods. Nevertheless, current experimental 
techniques have developed to the point that they are on the verge of 
reaching the region of these intrinsic limitations in earth-bound experiments. 
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4. Conclusions 

In this report we have shown that gravity causes experiments to 
become inaccurate when the critical point is approached sufficiently 
closely. We have considered the ranges in AT* and Ap* which are inaccessible 
in a variety of earthbound experiments and how these ranges shrink when 
gravity is reduced. Some of our conclusions are summarized numerically 
in Table 4.1. 

The critical point is a focal point of anomalous behavior of many 
physical properties. In order to understand the nature of these anomalies 
it is desirable to approach the critical point as closely as possible. 

In fact, the most recent experiments (Hocken and Holdover, 1976; Balzarlni 
and Ohrr, 1972) indicate that at temperatures as close a 1 mK from che 
critical temperature the anomalous behavior near the critical point of 
fluids, although appearing to approach Ising mode' behavior still 
differs from the behavior predicted by some theories for that model, k 
disappearance of this difference at the critical point would have 
considerable theoretical significance. 

Experiments In a low-g environment will provide opportunities to 
study the nature of the anomalies in a range of temperatures near the 
critical temperature inaccessible in earthbound experiments. To 
discuss which anomalies are best suited for study in ^pace, it is 
convenient to classify aromalles as "strong" and "weak” (Griffiths ard 
Wheeler, 1970). In this rough classification, quantities which behave 
like (T-T c ) ^ on the critical lsochore are called strongly anomalous if 
p is greater than, perhaps, 0.5. Thus the isothermal compressibility, 
the constant pressure specific heat, the thermal diffusivity, and the 
correlation length exhibit strong anomalies in pure fluids. In contrast, 
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properties for which p is less than 0.5 (or even negative) exhibit weak 
anomalies. Examples in pure fluids are: the velocity of sound, the constant 
volume specific heat, the dielectric constant, and the shear viscosity. 

In general, strongly anomalous quantities can be studied sufficiently 
well on earth (particularly with refinements of optical techniques 
mentioned In Chapter 2) so that low-g experiments do not seem desirable at 
present. If the technology advances to the point where the non local 
effects of Chapter 3 are encountered, this statement would have to be 
revised. On the other hand, those properties which exhibit weak anomalies 
and which must be studied in bulk samples are excellent candidates for 
scientifically valuable low-g experiments. The study of these weak 
anomalies In pure fluids is hampered by the strong variation of fluid 
density with height on earth. 

The weak anomaly In the index of refraction is a particularly 
important one for low-g studies. Thorough understanding of the index of 
refraction anomaly is of Interest in Itself; however, it is of even 
greater interest in Its Impact upon the interpretation of equation of 
state data obtained by optical techniques. Other good candidates for 
study in a low-g environment are the C^ and viscosity anomalies which, 
with present technology, must be studied in bulk samples at least several 
tenths of a millimeter high. 

In the areas of phase transition phenomena in fluids, low-g studies 
of spinodal decomposition, nucleatlon, and macroscopic phase separation 
all seem appropriate. At this time, it is hard to present a detailed 
justification of low-g studies because the phenomenology is not as well 
developed in these areas as It is for the thermophysical properties. We 
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believe exploratory experiments which exploit the long sedimentation 
times available in spacelab are appropriate. Ue agree with the Overstudy 
Committee (Dodge et al . , 1975) that very simple experiments will be 
extremely Important. The Overstudy Committee described one such experiment 
(Dodge et al . t p.9) which we heartily endorse: 

Take a set of samples of a one-component fluid with different 
densities. Included should be a sample at subcritlcal density, a 
sample at the critical density and a sample at supercritical 
density. Let the temperature vary either continuously or dis- 
continuously through the transition temperature and take photo- 
graphs at regular temperature intervals. Compare these photographs 
with those obtained at esrth under the same instrumental conditions. 
Conduct similar experiments for a set of samples of a binary mixture 
at concentrations smaller than, equal to and larger than the 
critical concentration. 
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Appendix A. Parametric equation* of state for fluids near the critical 
point . 

A.l Introduction 

Let p be the density. V the volume, P the pressure, T the temperature, 
w the chemical potential, ■ p ^(3p/3P)^ the isothermal compressibility 
and Cy the heat capacity at constant volume. It will also be convenient 
to introduce a "symmetrized" compressibility 

Xy * (t)p/3p)y s P*Ky (A.l) 

The thermodynamic properties are made dimensionless by expressing them 
in units of appropriate combinations of the critical temperature T c , the 
critical density p c and the critical pressure P c - Specifically, we 
define 

T* = T/T c , p* = p/p c , P* = P/P c 


(A. 2) 

= ^P c / p c » xf = X T p c /Pc* c v = C V T c /VP c 

Note that the reduced specific heat C* is taken per unit volume. In 
addition we define the differences 

An* = Mp.T) - w(p c ,T)}p c /P c 
AT* = (T-T c )/T c 
(p-p c )/p c 


Ap* = 
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(A. 3a) 
(A. 3b) 
(A. 3c) 



The behavior of various thermodynamic properties is described by 
power la 1 78 when the critical point is approached along specific paths in 
the AT* - Ap* plane. One customarily defines four critical exponents 
a, 8, y, and 6 associated with the thermodynamic behavior of .lulds near 
the critical point. (Rowlinson, 1969). The exponent a describes the 
divergence of Cy along the critical isochore. 

(T*T C , p = p c ) C*/T* * £i(aT (A. 4a) 

the exponent 6 characterizes the shape of the coexistence curve 

(T sT c ,p = P cxc ) Ap* xc = * 3 UT*| 6 , (A. 4b) 

the exponent y describes the divergence of the compressibility 
(T * T c , p = p c ) Xf s r (aT*P (A.4c) 

and the exponent 6 specifies the behavior of Ay along the critical 
isotherm 

(T ■ T c ) Ay* = D (Ap*)|Ap*f 6-1 (A.4d) 

The description of the thermodynamic behavior is based on a scaling 
hypothesis (Widom, 1965; Fisher, 1967; Levelt Sengers, 1974, 1975). This 
scaling hypothesis suggests that the equation of state upon approaching 
the critical point will be asymptotically of the form (Griffiths, 1967) 
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&p* |ap *| 6 ^h(x) 

( A.S) 

aT*/|ap*| 1/b 

(A. 6) 


The scaling hypothesis implies that the thermodynamic critical exponents 
introduced in (A* 4) satisfy the relations 


2 - a = 3(6+1 } 

y = 6(6-1) 


(A. 7) 


so that only two exponents can be chosen independently. 

The function h(x) in (A. 5) must satisfy a number of stability and 
analyticity conditions (Griffiths, 1967' . However, xt turns out to be 
very difficult to formulate an explicit mathematical expression for the 
function of h(x) that would satisfy all required analyticity conditions 
and which would be analytically integrable to calculate the specific 
heat (VicentJni-Missoni et al ., 1969; Schmidt, 1971; Levelt Sengers 
et al ., 1976) * 

These problems are solved by using parametric equations of state 
(Josephson, 1969; Schofield, 1969). In the parametric equations of 
state the relationship between the physical variables Ay*, AT* and Ap* 
is given implicitly via two parametric variables r and 6. The constraint 
that the scaling law (A. 5) is satisfied is met by the fcllowing choice 
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Ay* • r 66 H(e) 

aT* * r T(e) (A. 8) 

Ap* s r® M(e) 


The variable r is meant, in some sense, to describe a distance from the 
critical point and the azimuthal variable 0 a location on a contour of 
constant r. The cri^eal isochore, the critical isotherm and the 
coexistence curve are all curves of constant angle 3. In Fig. A1 we 
3hov tho location of these curves in the Ay* - aT* plane and indicate 
the meaning of the parametric variables r arid 6. 

Parametric equations of state with various choices for the functions 
H(e), T( e) and M( 0) have been used successfully to represent experimental 
equation of state data in the critical region of fluids. The parametric 
equations of state used in this repot are the restricted linear model, 
the restricted cubic model and the Wilcox-Estler equation of state. 
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Fig. Al. Coexistence curve, critical isochore and critical isotherm in 
the plane with coordinate axes Av and t = aT*. 
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A. 2 Restricted linear model. 


The restricted linear model equation of state corresponds to the 
choice (Schofield, 1969; Schofield, Lltster and Ho, 1969) 

dp* = a 1 r B 6 e(l-e 2 ) 

AT* = rd-b^e 2 ) 

Ap* = k-| 

with 

b 2 * ^ (A.9d) 

1 v 6 - 1 ) 0 - 28 ) 

where a^ and are adjustable constants. In this formulation 0 = 0 on 
the critical isochore. 9 * ± 1/bj on the critical isotherm and 9 = = 1 on 
the coexistence curve. The sign of 9 corresponds to the sign of £p*. 

The restricted linear model yields for the compressibility 

and for the singular contribution to the specific heat 
* 

Suiisa = a.k. r - (a. id 

T* 1 1 2a( 6-3) 

so that a contour of constant r may be interpreted as a contour of 
constant anomalous specific heat (Schofield, Lltster and Ho, 1969). 

The restricted linear model has been used to fit experimental data 
fcr magnets (Schofield, Lltster and Ho, 1969) and fluids (Hohenberg and 


(A. 9a) 
(A. 9b) 
(A. 9c) 
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Barmatz, 1972; Huang and Ho, 1973; Thoen and Garland, 1974; White and 
Maccabee, 1975; Levelt Sengers and Sengers, 1975; Levelt Sengers et al., 
1976). It can be Integrated and also be fitted tc the pressure (Murphy 
et al ., 1973, 1975). 

Restricted linear model parameters for a variety or fluids are 
presented in Appendix D. 

A. 3 Restricted cubic model . 

The restricted cubic model equation of state corresponds to the 
choice (Ho and Litster, 1970; Huang and Ho, 1973) 

An* = a 2 r 36 e(l-e 2 ) (A. 12a) 

at* = r(l-b 2 2 e 2 ) (A. 12b) 

Ao* = k 2 r 6 e(l+c e 2 ) (A. 12c) 

with 

b 2 = — , c = (A. 1 2d) 

C 3-28 3-2B 


where a 2 and are adjustable constants. Just as in the linear model, 
d • 0 on the critical isochore, 0 - 1 l/b 2 on the critical isotherm and 
0 * 1 1 on the coexistence curv^ , 

The restricted cubic model yields for the compressibility the 
simple form 



(A. 13) 


so that a contour of constant r corresponds to a contour of constant 
compressibility • 
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When Che restricted cubic aodel la fitted to experimental data the 
quality of the representation is comparable to that of the restricted 
linear model. Restricted cubic model parameters for a variety of fluids 
are presented In Appendix 0. 

A. A Wilcox-Estler model . 

The parametric equation Introduced by Wilcox and Estler (1971] has 
been usad to analyze density gradient profiles in the extreme vicinity 
of the critical point (Estler et al ., 1975; Hocken and Moldover, 1976). 
It is defined through the equations 

s„*= 2 r Y Y e (r,»)H(ei (AJ4a) 

Y 


iT* = r 0 


(A. 14b) 


with 


Ap* = Y B (r,e) [1 + B W(e)] 


W(e) 


1 - e/e > 

1 - #/ 0 , 


(A. 14c) 


(A. I4d) 


Y(r.e) = Y q r (l-e/e^' 


a = 1 - e 0 /e x 


(A. 14e) 
(A. 14f ) 


This equation of state is designed such that 


X = r ~ y 


(A. 15) 


so that a contour of constant r - (1/x^. ) ^^corresponds to a contour of 
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constant compressibility as in the restricted cubic model. It ha 9 the 

additional feature that AT* is a linear function of both r and 0 . 

In this formulation lines of constant 6 are straight lines as Indicated 

in Fig. A2. In this model 6 - o corresponds to the critical Isotherm , while 

9*0 on the critical Isochore and 6*0 on the coexistence boundary, 
o * 

For the special choice 0A “ 3/2 the Wilcox— Estler model reduces to the 
restricted cubic model (Estler etal. , 1975). 
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CRITICAL POINT 


Fig. A2. Coexistence curv°, critical isochore and critical isotherms in the 
plane with coordinate axis r = (1/xy )^ Y and t = aT*. 
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Appendix B. Calculation of density profile . 

We assume that at each level z the local chemical potential 
u(p(z),T) equals tne chemical potential of a system with uniform density 
p - p (z) at temperature T in the absence of gravity. The conditions 
under which this assumption is valid are discussed in Section 3 of this 
report. Since in the presence of gravity the total chemical potential 
is the sum of the local chemical potential v(p(z),T) and the gravitational 
potential mgz, we have 



u(d(z),T) - u(o(z 0 ),T) = - mgAz (B.l) 


where m is the molecular mass, if the chemical potential is taken per 

particle and where Az = z - z . For convenience we take the reference 

o 

level z as the level where p « P c * In terras of dimensionless quantities 
we write 


Az = - H Ay* 


(B.2) 


with 


H - P c /(p c mg) (B.3) 

The quantity Ap* is defined in (A. 3a) and the product p^ m is the mass 
density at the critical point. The quantity H represents a scale height 
for the chemical potential of a fluid in a field of gravity g. H q 
represents the scale height for the chemical potential in the earth’s 
gravitational field where g - g . In Appendix D we present the scale 
factors for the critical Isotherm for a number of fluids in the earth’s 
gravitational field* where 
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(B . 4 ) 


D H 

o 


DP c /p c ® 8 


o 


(These scale factors for the critical isotherm differ from H q by a 
numerical factor, D, which varies between 1,2 and 3.2 among various 
fluids). Introducing a gravity ratio 


g* ** g/g Q 


(B.5) 


we may rewrite (B.2) as 


AZ 



(B.6) 


Substitution of the scaled equation of state (A. 5) into (B.6) 
yields the relationship between height and density at a given temperature. 
In practice we either use the restricted linear model (A. 9), the restricted 
cubic model (A. 12), or the Wilcox-Estler model (A. 14). 

It follows from (A.l) and (B.6) that the density gradient 
(ap*/3z) T is given by 



(B.7) 


In practice the density can be determined as a function of height 
using a float densimeter (Greer et al . , 1974) or by measun. ~ the capacitances 
between a stack of horizontal conduction plates (Weber, 1970). The 
density resolution in such experiments is limited by the variation of 
the density over the height h of the float or the distance h between the 
two capacitor plates. If the local density is to be obtained with a 
precision p we must require that 
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lp*Uz+h/2) - p*(&z-h/2)l a h|f§*! s P 


(B.8) 


It thus follows from(B.7) that errors due to gravity in determining 
densities are avoided if 

* H 0 p 

X T S g*!? (B.9) 

or in terms of the cubic model parameters, using (A. 13), 


, »/sa *\ ,/Y 

\H oP a 2 ) 


For xenon this condition reads 


r s 6.2x10' 4 ($!!l) 0 * 84 


(B. 10) 


(B.n ) 


where the height h is to be expressed in m. In particular at the critical 
isochore p - p 

c 


AT* 2 6.2 x 10" 4 (3^) 0 - 84 


(B. 12) 


and at the critical isotherm T ■ T 

c 


Ap*i 0.076 (2^1) 0 - 3 


(B.13) 


76 



1 


t 




Appendix C. Calculation of correlation length . 

C.l. Definition of correlation lenp/ch . 

When a fluid approaches the gas-liquid critical point, lt6 thermodynamic 
state is accompanied by large fluctuations in the density. The magnitude 
and spatial character of these fluctuations are described in terms of a 
correlation function defined as (Stanley, 1971) 

P 2 G(|R-R'|) = < p(R)p(R’) > - P 2 (C.l) 

where p (R) is the local (number) density at position $ and 
p the average equilibrium density which is independent of the position 
^ (not considering the presence of external forces such as gravity). 

The zeroth order moment of the correlation function is related to the 
isothermal compressibility by the fluctuation theorem 

k R T x T - P 2 /dR G(R) (C.2) 

The correlation function G(R) ■ G(Ap*,AT*;R) is a function 
Ap* and AT* as well as of R (where R ■ jR - R * | ). 

The spatial extent of the fluctuations is characterized by a correlation 
length C- It is defined as (Fisher, 1964, 1967) 

£ - 2 _ 1 fdR R 2 G(R) 

' " 6 ~ (C. i) 

JdR G(R) 

The correlation length § diverges at the critical point. In particular 
along the critical isochore it follows the power law 

(T i T c , p - p c ) S - f 0 (AT*)' V (C.4) 
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Another exponent n is introduced to specify the nature of the 
dependence of the correlation function on the distance R. It is defined 
such that at the critical point Ap* = 0, AT* * 0 


G(o,o;R) = 


,1+n 


The exponent is zero in the classical theory of Ornstein and Zernike 
(Fisher, 1964). In current theories n is small, but finite. 

The correlation function exponents v and n are related to the 
thermodynamic critical exponents a, 6, y,6 introduced in Appendix A by 
the relations (Fisher, 1967; Widom, 1974) 


Y = -(2-n) 

(C.6) 

3v = 2-a 

(C.7) 

2-n _ 6-1 

T" 6+T 

(C.8) 


The* relations (C.7) and (C.8) are sometimes referred to as hyperscaling 
relations (Levelt Sengers and Sengers, 1975). 
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C.2 Correlation length as a function of density and temperature . 

The correlation length £ ■ £(Ap*»AT*) is a function of density and 
temperature. The hypothesis of scaling for the thermodynamic behavior 
can be extended to the correlation function G(Ap*» Al*;R) as a function 
of Ap*. AT* and R (Kadanoff, 1966). The scaling hypothesis implies that 
the correlation length £ can be written in the form (Sengers and Levelt 
Sengers, 1976) 

1 

5 s S 0 R(x/x Q ) (r _1 x T *) n (C.9) 

where x * AT*/ | Ap* l^^is the thermodynamic scaling variable introduced 
in (A. 6) and x q - - AT*/ ) A P* XC I B * , where B is the amplitude of 

the power law (A. 4b) for the coexistence curve. The function R(x/x o ) is 

St 

a universal function such that R(“) * 1 at the critical isochore x * *. 
The amplitude n be deduced from light scattering data (Chu, 1972). 

In this report we have calculated the correlation length in the 
approximation that R(x/x q ) * 1 independent of x 

1 

{ • 5 0 (r _1 x T *) ^ (C.10) 

This approximation, though not strictly valid, is adequate for the 
purpose of this report (Sengex and Levelt Sengers, 1976). In this 
approximation a contour of constant compressibility y* coincides with a 
contour of constant correlation length ?. For this reason we found it 
convenient to calculate £ using the restricted cubic model equation of 
state (A. 12) for which 
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V 


(C.lla) 


C = 


r" 


AT* -■ r(l-b| e 2 ) 


a p* ■ kg r B 0(1+ c e 2 ) 


with b 2 and c again given by (A.12d). 


(C. 1 Id) 
(C.llc) 
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C.3 Correlation length as a function cf height . 


Once the density profile and density gradient profile has been 
calculated from (B.fc) and (B.7), the correlation length as a function of 
height follows f^tu, (C . 10) . In terms of the parameters of the restricted 
cubic model the relationship between the correlation length £ and the 
height 2 at a given temperature oT* is determined by 

5 = s„ UT*r P - bj «T (C.12.) 

, jiT*| 86 P - b| e 2 r 66 9(l-e 2 ) (C.12b) 

9 

The correlation length assumes its max-uaum value at the level 
Az * 0 (6 = 0) corresponding to the critical density. The height at 
which the correlation length will be reduced by a factor 1-p is determined 
by the conditions 





o-p; 






(C. 1 3) 


so that 


Az 


. 5. S A* ^ V ’ 


2 2 
v b 2 - b 2 B6p 


) 


(C.1A) 


H 0 a ? p : a 

: *7i r (-)*h T *i‘ ,w 

9 D 2 ' 

The rate at wnich the correlation length f , at arbitrary le • 
vary as a function of height is Jetermined by 
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(C.17b) 


T he correlation length is the fundamental length scale which 
determines the anomalous behavior of thermodynamic ard transpot t properties 
near the critical point. In actual experiments these properties are- 
measured as averages over a finite height h. Such experiments will be 
reliable when the correlation lenr^h over this height does not vary 
within the precision p desired. Thus in order to avoid errors due to 
gravity we require that 


t 
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|e(Az + h/ 2 ) - (az - h/ 2 )| s p cUz) (c. 18 ) 


which ve approximate as 


1 / d£ \ p 

5 V 3z /r s " 


( C . "s 9 ) 


Upon substituting (C.lls) and (C.16) into (C.19) we conclude that the 
error.* 1 due to gravity will become appreciable unless 


r"" , ) 


2 

2 b 2 ? 


V 


266 b 2 2 8 v!- n2 ) + (^-3e : )(l-b 2 2 r : ) 


(C.20) 


When combined with (C.llb) and (C.llc) this equation defines a range in 
the AT* - Ad* plane where the measurements become inaccurate due to 
gravity effects. In particular at the critical temperature AT* » 0, 
this condition becomes 


where D is the amplitude of 
t -rticular for xenon 


.1/: 


(■* 0 ° “) 


iC . 21 ) 


the critical Isotherm defined in (A.4d). In 
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V 


(C.22) 


!ip*: i 0.16 l 2 *-) 0 * 23 
u 


where h Is to be expressed in meters. 

At the critical isochore Ap* * 0 at T + T c> (H/3*)^ • 0 and (C.19) 
is not a good approximrtion to (C.18). The distance Az^ over which the 
correlation length i9 reduced ty a factor 1-p is given by (C.14). In 
ord~r to avoid gravity effects we must require that this distance be 
larger than the height (strictly h/2) over which the properties are 
measured 



(C.23) 


It thus follows from (C.l4> and (C.23) that at the critical isochore 


kT*< * fen *1 { v } 1/2 


L 2H o 


P 


l/c-5 


(C.24) 


In particular for xenon 


j AT*j. 1.85x1 O' 3 /3^L\°- 65 
VP ! 


(C.25) 


where h is to be expressed in meters. 
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Table C.l. Critical Ration Parameters for a Nuabcr of 


Fluids AasuBlna Universal Effective 


Restricted 
Cubic Model 


Critical Point Parameters I 


P, 

p.a 



c 

MPa 

C 3 

kq/m 3 


*He 

0.11678 

41.45 

-3.3099 

4 He 

0.22742 

69.6 

5.1895 

Ar 

4.865 

535 

150.725 

Kr 

5.4931 

908 

209.286 

Xe 

5.8400 

1110 

289.734 

p-h 2 

1.285 

31.39 

32.935 

«2 

3.398 

313.9 

126.24 

°2 

5.043 

436.2 

154.580 

h 2 o 

22.06 

322.2 

647.13 

020 

21.66 

357 

643.89 

co 2 

7.3753 

467.8 

304.127 

nh 3 

11.303 

235 

405.4 

SF 6 

3.7605 

730 

318.687 

CHj 

4.595 

162.7 

190.555 

c 2 H 4 

5.0390 

215 

282.344 

C2 H 6 

4.8718 

206.5 

36 .33 

C 3 H 8 

4.247 

221 

369.82 

Notes : 

: o*0. 

100 S - 

0.355 


y m 1. 

190 « - 

4.352 


v * 0.633 n * 

0.12 


Restricted 
Linear Model 

k l 

*1 

0.924 

4.58 

1.021 

6.40 

1.309 

16.1 

1.309 

16.1 

1.309 

16.1 

1.156 

9.6 

1.361 

18.2 

1.309 

15.6 

1.622 

21.6 

1.622 

21.6 

1.436 

21.3 

1.573 

21.4 

1.337 

23.9 

1.361 

17.0 

1.350 

17.5 

1.416 

20.2 

1.451 

20.2 

b l"- 

1.3909 


4 •? 


0.818 4.05 
C.904 5.66 

1.160 14.2 
1.16C 14.2 

1.160 14.2 
1.024 8.5 
1.206 16.1 

1.160 13.9 
1.438 19.1 
1.438 19.1 
1.273 18.9 
1.394 19.1 
1.185 21.2 
1.206 15.1 
1.197 15.5 
1.255 17.9 
1.286 17.9 


b 2 “ * 1.3100 


906 

946 

2260 

1500 

1310 

10410 

2560 

2790 

8980 

7950 

3470 

7134 

1730 

6260 

5540 

5520 

3830 


The lest quantity tabulated, DH » D P /to. ng_) il the scale factor f*»r the density 

O £ *• *■ 

vs height profile at the critical teaperature. D la found from the relation: 

* 1 2 2 

D • a.k. b. (b, - 1)» which is taken from Sengers et ale [1976]. 



















Appendix D. Physical constants for various fluids . 


In this Appendix we present critical parameters and critical region 
parameters Tor a number of fluids* taken from r survey prepared by 
Sengers and Levelt Sengers [1977]. The restricted linear model and 
cubic model parameters in this table represent Informed estimates assuming 
universal effective critical experiments corresponding to the range 



Appendix E. Parameters used for xenon in this report 


Critical point parameters .(Levelt Sengers et al., 1976) 

P - 5.8400 MPa 
c 

mp c - 1110 kg/m 3 

T « 289.734 K 
c 

H = P /p g * 536 m 

O C CO 


Restricted linear model (Levelt Sengers et al., 1976) 


a - 0.100 

A - 1.71 

k x - 1.309 

8 = 0.355 

B = 1.827 

a 1 = 16.1 

Y = 1.19 

r * 0.0813 

bj 2 = 1.3909 

6 = 4.352 

D * 2.44 


Restricted cubic mode" (Sengers and Levelt Sengers, 1976) 


a = 0.100 

A « 1.68 

k 2 - 1.160 

3 = 0.355 

B - 1.827 

a 2 - 14.2 

Y = 1.19 

T = 0.0817 

b 7 2 * 1.3100 

6 = 4.352 

D « 2.44 

c * 0.0393 

Vilcox-Estler model (Hocken 

and Moldover, 1976) 


n = 1.1379 

8 - 0.3293 

0 = 0.1076 

c 

o 

n '= 0.14 

Y - 1.229 

Y d = 0.323 

c 

A * 4.4 

o 


Correlation length (Swinney and Henry, 1973) 

v = 0.633 £ - 2.0x10 

o 
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a r hT 


Nomenclature 


A 

as 

amplitude 

of 

power law for specific heat C v 

8 1 

- 

parameter 

of 

linear model 

a 2 

- 

parameter 

of 

cubic model 

B 

tz 

amplitude 

of 

power law for coexistence curve 

b l 

= 

parameter 

of 

linear model 

b 2 

= 

parameter 

of 

cubic model 

C 

V 

e 

constant-volume heat capacity 

C * 

V 

- 

C T /V P 
v c c c 




c * parameter of cubic model 

D - amplitude of power law for chemical potential along critical 

isotherm 

d ■ sample thickness 

g(R) c correlation function as a function of distance 

g « gravitational acceleration 

2 

g^ = 9.80 m/s - earth's gravitational acceleration (units) 

g* *= g/g Q « gravitation acceleration relative to its value at the 

earth’s surface 
h * sample height 

H = P c /p c mg * scale height for chemical potential 

H q * P c^ p c m ®o r scale height for chemical potential on earth 

k “ wave ve . r 

p -1 (ap/aP) T - isothermal compressibility 
Boltzmann's constant 
k^ * parameter of linear model 

k 2 = parameter of cubic model 

m * molecular mass 

n ** refractive index 
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p 

p 

c 

p* 

AP* = 
P 

r = 

-+ 

R 

R 

R 

S 

T 

T 

c 

T * 

AT* = 

t 

V 


z 

z 

o 

Az 


refractive index at critical point 
(on/3p*)^ 

* value of at critical point 
Pressure 

critical pressure 
P/P 

c 

(p-p)/p 
c c 

desired precision 
parametric variable 

spatial coordinate 

magnitude of spatial coordinate 

scaling function for correlation length (Appendix C.2) 

structure factor 

temperature in Kelvin 

critical temperature 

T/T 

c 

(T-T )/T 
c c 

AT* 

volume 

1/8 

scaled thermodynamic variable AT*/jAp*| d or scaled wave vector 

b-i/b 

parameter of Wilcox-Estler model 
height 

height at which p * p c 
z-z 
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F 

\ 


Greek letters 


exponent of power la\' for specific heat C y 
exponent of power law for coexistence curve 
exponent of power law for compressibility 
amplitude of power law for compressibility x T 
decay rate of entropy fluctuations 

exponent of power law for chemical potential along critical 
isotherm 

dielectric constant 

3 / r - parameter cf Wilcox-Estler model 
ox 

exponent for spatial dependence of correlation function 
at the critical point 
parametric variable 

value of 0 on critical isochore in Wilcox-Estler model 
value of 0 on coexistence curve in Wilcox-Estler model 
correlation length 

amplitude of power law for correlation length 
wave length 

chemical potential per particle 

up /P 
c c 

{ W (p >T) - u(P c ,T)}P c /P c 

exponent of power law for correlation length 
number density 


critical density (particles/volurae) 
average density 


(p - p c )/p c 
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f 


t 

> 


cxc 

* 

) 

cxc 


density of vapor or liquid at coexistence 


P /p 
cxc c 


cxc 


(p„ - p„)/p„ 

cxc c c 


cross section area 


turbidity 

(9p/3|j) T = p 2 K t[1 

(3p*/3y*) « x T p c /P c " 


frequency 


i 


* 
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